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We consider the renormalization of the bending and Gaussian rigidity of model membranes induced 
by long-range interactions between the components making up the membrane. In particular we 
analyze the effect of a finite membrane thickness on the renormalization of the bending and Gaussian 
rigidity by long-range interactions. Particular attention is paid to the case where the interactions 
are of a van der Waals type. 
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I. INTRODUCTION 



In the Monge gauge the simplest Hclfrich Hamiltonian used to describe the bending energy of a membrane, 
whose height fluctuations in the direction perpendicular to the projection plane are written as z = /i(x), is given by 

H Hel = ^ f dx(V 2 Mx)) 2 . (1) 



Here Kb is the bending rigidity and it is principally generated by short range lipid- lipid interactions. Clearly however Kb 
will also depend on long-range interactions between the membrane components such as van der Waals and electrostatic 
interactions. Unlike short range interactions, these long-range interactions can in some cases be modified rather easily, 
for example one can screen electrostatic interactions by adding electrolyte to the system. 

A number of authors have considered the problem of renormalization of bending and Gaussian rigidities due to long- 
range interactions @, 0, & 0> El 13 ■ Clearly this renormalization is important in determining membrane morphology 
and is thus of great physical and biological importance. The first effects studied were those due to the presence of a 
surface charge on the membrane surface, more precisely an inner and outer surface charge. These systems were first 
studied using the Poisson Boltzmann mean field theory in its linearized form |2( and subsequently in its full non- linear 
form In the case of fixed symmetric surface charge (with respect to the two surfaces of the membrane) it is found 

in these studies that renormalization of the bending rigidity is positive, however the renormalization of the Gaussian 
rigidity is found to be negative. At the mean field level, even in the presence of a screening electrolyte solution, the 
effective interaction between like charged surfaces is repulsive, the increase of the bending rigidity is thus easy to 
understand from a physical point of view. The reduction of the Gaussian rigidity in these systems is less obvious to 
understand from a more intuitive stand point. The general conclusion coming from the studies of charged membranes 
is that the presence of surface charges can favor their buckling and perhaps induce the formation of spherical vesicles. 

The literature on the renormalization of membrane rigidity by long-range interactions contains a few contradictory 
results, some of which this paper will try to partially resolve. Firstly a number of authors predict that non retarded 
van der Waals interactions between the membrane components will lead to a reduction of Kb 0,0- These non-retarded 
interactions are attractive and the general formalism developed by Netz |7j predicts that attractive pair-wise inter- 
membrane interactions will lead to a reduction of Kb- This seems eminently reasonable from a physical perspective, 
attractive interactions favor points on the membrane becoming closer together and thus a flat sheet which is rolled 
into a cylinder should have a lower energy. This reasoning is however rather deceptive, indeed we shall later see that 
while attractive forces favor the formation of a cylinder from a plane sheet, they may actually inhibit the formation 
of a spherical vesicle. 

Taking a flat membrane and rolling it up to form a cylindrical surface of length L and radius R costs a bending 
free energy given by 

H h = (2) 

As mentioned above, it seems clear that at a fixed surface area, if the interaction is attractive then it is energetically 
favorable to form a cylinder from a flat sheet if the only energies involved are due to attractive long-range interactions 



between the membrane components. Hence, if this is the case, the renormalization of Kb above should be negative. 
In a recent paper Q we showed that this is indeed the case, in agreement with Netz who used his approach for a 
generic pairwise potential applied to the case of effective dipole-dipole interactions. In [3] the static van der Waals 
interactions are taken into account via the difference in the dielectric constant of the membrane of finite thickness 
and the surrounding bulk, e.g. water, using the Lifshitz theory 9]. If S is the membrane thickness, we showed that 
for large R and when the dielectric constant of the membrane cm is close to the dielectric constant of water ew, the 
van der Waals contribution to the free energy can indeed be seen as a renormalization of the bending rigidity i.e. the 
free energy is given by Eq. J2J with 

K b = - — ^[31og(A5) + 0.02954], (3) 

where A = (ew — £m)/(cw + £ m) and A ~ 1/a where a is a microscopic cut-off scale. In an older paper |4| the same 
system was analyzed and the resulting Kb was predicted to be positive, implying a membrane stiffening due to the 
van der Waals interactions. More precisely in Q| it is predicted that, for large R, Hb ~ Const, x hx(RA)/R. The 
difference between the results of and can be traced to the fact that the expression for the free energy given in 
does not go to zero when the thickness S of the membrane goes to zero, i.e. when no membrane is present. When 
the result of is normalized to give zero when 6 = then the result of is obtained. The apparent logarithmic 
behavior in the result of £| is found to disappear, although the mechanism is very subtle and involves some very 
lengthy analysis. The calculation of § also predicts that the renormalization of Kb is negative for arbitrarily large 
differences in dielectric constants, i.e. beyond the dilute approximation where the description in terms of an effective 
pairwise interaction is no longer valid. 

An electro-neutral membrane of infinitesimal thickness containing monovalent charges, such as a 1-1 electrolyte, 
has been studied in [(| . In this study the surrounding medium is taken to be a non-ionic solvent (a somewhat unlikely 
situation) . It was shown that for this system a cylindrical geometry has a bending free energy due to the charges of 

H b = -^\n(R/\ D ), (4) 

here 1/Xd = StthqIb where Ib is the Bjerrum length and no is the (positive/negative) surface charge density of the 
membrane. This predicts an R dependent bending rigidity 

n b = -^HR/\ D ). (5) 

Subsequently Netz analyzed this problem in his general formalism where the area elements of the membrane interact 
via a generic pairwise interaction. In order to do this he was obliged to calculate and effective pairwise interaction 
for the system. He found that for this particular system 

Kb = ~384^ (6) 

It is interesting to ask what is the origin of these differences. In the method of [f| one looks at the fluctuations 
about the mean field and effectively calculates a functional determinant related to the surface charge fluctuations. 
In this calculation the energy of a flat membrane is explicitly subtracted off to yield the bending energy. In the 
Netz formalism the energy is expanded to quadratic order in small height fluctuations about a flat membrane, the 
perturbation being taken about the true area of the membrane as opposed to the projected area. Possible reasons for 
the observed discrepancy, assuming both calculations to be formally correct, are 

• There is an ensemble difference: the energy of a flat membrane is only explicitly subtracted off in the calculation 
of@ 

• As pointed out by Netz 0, his calculation is done as a perturbation about a flat membrane but the calculation 
of [6j is explicitly done for a cylindrical (and spherical geometry). This could lead to a difference. 

• The effective interaction needs to be calculated in this system and there could be a difference in the approxi- 
mations used. 

In what follows we shall show that the formalism of Netz can be applied with an explicit subtraction of a flat 
membrane energy of same area and that, under rather weak assumptions, the same result is recovered. This cannot 
therefore be at the origin of the disagreement. In addition we show that for a general potential V between membrane 
components a calculation done in a cylindrical geometry agrees with the result of Netz based on a perturbation about 



2 



a flat plane, however the potential V must be of sufficiently short range. This result makes sense as the low momentum 
height fluctuations correspond to regions of low and slowly varying curvature, and it is thus unlikely that this could 
lead to a difference, at least for short range potentials. We then use the geometric method of calculation to compute 
the Gaussian rigidity K g by examining the case of a sphere. In some cases the formulas for Kf, and K g obtained for 
an infinitcsimally thin membrane need to be regularized by a short distance cut-off in the potential V. However 
membranes have a finite thickness and it is thus interesting to investigate the effect of a finite membrane thickness 
to see if this regularizes the rigidity renormalizations. The renormalization of rigidities for finite membrane thickness 
are therefore also derived via the geometric route. 

As an application we carry out the calculation of renormalization the bending and Gaussian rigidities of a membrane 
of finite thickness S having a differing dielectric constant to that of its surrounding medium. . This calculation is carried 
out in the dilute limit where the two dielectric constants are close to each other and is equivalent to computing the 
effect of the zero frequency van der Waals interactions. In this calculation the effective pairwise van der Waals 
interaction is regularized with a real space cut-off, as opposed to the Fourier space cut-off employed in the eigen-modc 
expansion of . The result of this calculation with a real space cut-off is strictly identical to that of Eq. © up to 
the definition of the cut-off scale. The contribution to the Gaussian rigidity is also computed and found to be positive 
and with leading divergence in the small cut-off, again denoted by a, of the form S 2 /a 2 . 



II. THE NETZ APPROACH EXPLICITLY RENORMALIZED WITH RESPECT TO A FLAT 

MEMBRANE 

Here we revisit the approach of Netz for the computation of the renormalization of the bending rigidity of a close 
to planar membrane by long-range interactions. The argument we give is slightly different in that we use a formalism 
which calculates the energy of a fluctuating membrane with respect to that of a flat membrane with no fluctuations 
but with the same surface area. Reassuringly this approach gives exactly the same result. We imagine a membrane 
with projected area Aq in the Monge gauge over the region x E A a and of height h(x) in the z direction at that point. 
If the membrane is allowed to fluctuate in the z direction then the bending energy due to these fluctuations is given 
by: 



H b = I 
2 



AoxAg 



dxdx 1 y] I + (V/i(x))y I + (V/l(xO)V (V(x-x') 2 + (Mx)-Mx')) 2 ) 

- i^dxdx'V(y(x-x') 2 ) (7) 

Note that the dimensions of V as defined above are thus [23] /[Z 4 ], E denoting energy and and L length. In the case 
where V is a purely energetic interaction then is a purely energetic term However if V is an effective interaction 
induced by thermodynamic effects, (e.g. presence of electrolyte when the membrane is charged or static van der Waals 
interactions), then it will have a temperature dependence and Hb is thus strictly speaking a free energy. The real area 
(as opposed to the projected area) of the membrane is given by 

A= f dx Jl + (Vh{x)f (8) 
Ja 

and the second term in Eq. J7J) corresponds to the subtraction of the energy of a flat membrane of the same area, 
which can be thought of as a flat bulk membrane from which the membrane we study is drawn from. To quadratic 
order in the height fluctuations h we obtain 

H = - [ dxdx'VQx - x'|) - - [ dxdx'V(\x-x'\) 

2 JAqxAq 2 J AxA 

+ \J dxdx' V(\x - x'|) (VMx)) 2 + \j dxdx' V '^-j' l \ h( X ) - Mx')) 2 , (9) 

2 JA xAo ^Ja xA |x-x'| 

the second term depends on h through the area A which to this quadratic order is given by 

A = A + l [ dx(Vh(x)) 2 . (10) 

1 JAo 

This expression may be simplified in the limit of large A if the following integrals over all R 2 converge 

v = J dx V(\x\) = 2ir drrV(r) (11) 
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«i = / rix ^ff = 2tt / drV"(r) = 2tt(V(oo) - V(0)). (12) 



The above integrals converge when V(0) is finite, which can be ensured via a suitable regularization scheme, and 
when V(r) tends to zero quicker than 1/r 2 for large r. In this case we find 

H=\j dxdx'/i(x)G(x-x>(x'), (13) 

where the operator G is given by 

G(x) = -^V 2 ( 5(x)+^(x)-^M. (14) 

The two dimensional Fourier transform of G is given by 

G{q) = vi + ^!f- - 2n [°° drV'(r)Mqr) (15) 
* Jo 

where 



1 f 2 " 

Jo(q r ) = 7T" / d9exp (iqrcos (0)) (16) 
2tt Jo 

is a Bessel function of the first kind [l(J • This results may be written as 

G(q) = 2nj°°dr (l - £f - J (qr)\ V'(r), (17) 

recovering the result of Netz Q by using the identities 

2tt / drV'(r) = vi (18) 
Jo 

and 

/•OC 

7T / drrV(r) = -u 0) (19) 
Jo 

where the second assumption is valid providing that linv^oo V(r)r 2 — 0. An integration by parts can be used to 
rewrite Eq. as 



G{q) = 2vr / dr (q 2 r - 2qJ 1 (qr)) V(r). (20) 
Jo 

The small q expansion of Eq. (|20|l gives 

4 roo 

G(q) = ^- jf drF(r)r 3 + 0( g 6 ). (21) 
Going back to real space this low momentum term gives a contribution to the Helfrich energy of 



2 

where 



A Hel = S -^ f dx(V 2 Mx)) 2 (22) 

AoxAo 



8 Jo 

is the renormalization of the membrane bending rigidity. 



rOO 

8K b = o / <W(r)r 3 , (23) 
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III. GENERAL CALCULATION FOR CYLINDERS AND SPHERES 



In a cylindrical geometry the energy for an interaction V between membrane components with the corresponding 
energy of a flat membrane of the same area subtracted off is the bending energy. For sufficiently short range interactions 
it is given by 



H h = 



2irRL 



dz 



Rd6V 



'4i? 2 sin 2 (|) 



z 2 - 2tt 



dr rV(r) 



(24) 



The fact that the potential is short ranged is used in setting the limits of the z integration in the first integral at ±00 
and the limits of the r integration in the second integral (for the flat membrane) at and 00. In the first integral the 
only contribution that will be present for a short range potential will be from the region 9 = 0, more specifically the 
region where 9 = a/R and where a ~ 0(1). Making this change of variables and expanding the sin we obtain the 
first integral in the above to be 



To leading order in l/R this is 
For large R this can be written as 



dz 



wR 

dz I da 

ttR 



Till 



daV 



-irR 




R 2 



+ 0(l/R 4 )+z 2 



V{^a 2 + z 2 ) - 



2AR 2 Va 2 + z- 



- ; V'{\/a 2 + z 2 ) 



I = 



dr 



2tt 



d9 



-V(r) 



l («0 



24i? 2 



V'(r) 



where we have written a = rcos (9) and z = rsin(#). We thus find that to leading order 



with 



2tt 



R 



(25) 



(26) 



(27) 



(28) 



d9cos 4 {9) I dr r 4 V'(r) 
dr r 3 V(r). 



(29) 



Therefore we see that the calculation in an explicitly cylindrical geometry agrees with that obtained by the formalism 
of Netz. This makes it unlikely that the disagreement between Q and |y] is due to a difference of topology, at least 
for suitably short range potentials. However the long-range nature of the potential may play a crucial role in the case 
of a salty membrane with no external electrolyte. In the calculations presented above one sees clearly that the first 
order contribution to the bending energy comes from local curvature. This means that the energy associated with a 
given point comes from the local curvature about that point. For a sphere, for example, one is clearly neglecting the 
contributions from the sphere which come from the interactions with points on the opposite side. If the potential V 
has no intrinsic scale and has a power law form V(r) ~ l/r a at large r, then the contribution from a point due to 
points opposite to it is of the order R 2 ~ a where the term R 2 is the area term. This means that this contribution is 
negligible for a > 2, but plays an role for a < 2. Indeed one sees that the first order term in the energy of the system 
per unit area, which corresponds to the flat membrane is given by 2ir rdrV(r) and this is clearly divergent when 
a < 2. Furthermore it is clear from Eq. (|I2|) that even the analysis about the flat membrane breaks down in the case 
a < 2. 

Let us now consider the same calculation but for a sphere. For suitably short range potentials the bending energy 
of the sphere is given by 



H h = 2ttR z 



R' 



sm{9)d9 I d<f>V [R\/2 - 2cos (0)) - 2?r / dr rV(r) 
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An 2 R 2 



21! 



dr rV(r) — I dr rV(r) 



-4tt 2 R 2 / dr rV{r), 

J2R 



(30) 



where to obtain the above we have made the change of variables 2 — 2cos (9) = r 2 /R 2 . We see from Eq. <|3U[1 that if 
the potential V decays exponentially, or faster than exponentially, as a function of R, then the bending energy has 
an exponential, or faster than exponential, decay. For a sphere the bending energy takes the form 



Ht = 87TK;, + At: Kg, 



(31) 



where Kb is the bending rigidity and K g is the Gaussian rigidity • For short range potentials the above calculation 
implies that the presence of a long, but finite range, interaction does not renormalize the bending energy of a sphere 
to first order, rather we conclude that 5k 9 = —28Kb ~ so the renormalization of the Gaussian rigidity is of the opposite 
sign to that of the bending rigidity. We note that the signs of these results are in agreement with the results on 
charged membranes 00. It is interesting to note the calculations based on perturbations about a flat plane do not 
indicate any renormalization of the Gaussian rigidity by long-range interactions. This must be due to the inherently 
topological nature of the Gaussian bending energy which is insensitive to the geometry of the system and only 
depends on the membrane topology. It is thus normal that this energy is not picked up by a local perturbative 
analysis. Notice that if V is negative, and decaying at R — > oo sufficiently quickly, then the bending energy is actually 
positive - at variance with the intuition that attractive interactions favor bending ! 



IV. MEMBRANES OF FINITE THICKNESS 



In reality the membrane will always have a finite thickness. In the case of a cylinder we will examine the effect of 
long-range interactions on the bending rigidity of a cylindrical shell C of uniform thickness 5. We will take the outer 
radius of the cylinder to be at R + 5/2 and the internal radius at R + 5/2. If the cylinder is of length L then the 
volume of the shell is 2ttRL5. The energy of the cylinder is now given by 



H= \ [ dndr2V(|ri-r 2 | 

1 JCxC 



(32) 



Note that the potential V as defined above is now due to the interaction of volumes and not areas, it thus has physical 
dimensions [J3]/[Zi ]. The energy above can be written as 



H = itL J r\r?,dr\dr<idBdzV (^/(n - r 2 ) 2 + 2rir 2 (l - cos (0)) 



+ z- 



(33) 



where in the above the integration ranges are (R— 5/2, R + 5/2) for r\ and r 2 , (— 7T, tt) for 9 and (— oo, oo) for z. One 
now makes the change of variables 9 = a/ \jr\r-i, and recalling that r\ and r 2 are of order R, expand the argument of 
V to fourth order in a to obtain 



H = nL J y 1 r\T2dr\dr2dzda 



V (VOi -r 2 ) 2 + z 2 + a 2 ) 



a 4 V [y{rx-r 2 ) 2 + z 2 + a< 
24rir 2 ^(r x - r 2 ) 2 + z 2 + a 2 



(34) 



Again taking the limits of the a integrations to ±oo, carrying out the 9 integration, then replacing the coordinates 
(a, z) by the radial coordinates (r, 9') and writing = R + Xi for i = 1 and 2 we find that to order 1/R 

1 



H = 7r LR I dx\dx2rdrd9' 

V (^{x 1 -x 2 ) 2 +r 2 



(1 + 2R~( Xl+ X2) + W 2 * 1 * 2 ~ SR 2 ^ 1 + X2 ^ V 



X2j 



24R 2 



r 3 cos 4 ( 



\/{xi - xi) 2 + r 2 

Now carrying out the 9' integration and an integration by parts on the last term above we obtain. 



(35) 



H = Hi 



bulk 



H h 



(36) 
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where Hbuik is the bulk energy dependent only on the volume and given by 

H bulk = 2n 2 LR J d Xl dx 2 rdr V (y/( Xl - x 2 ) 2 + r 2 ) , (37) 



and Hf, is the bending energy given by 

H -** L 



d Xl dx 2 rdr (r 2 - 2(xi - x 2 ) 2 ) V ( y/fa - x 2 ) 2 + r 2 ) . (38) 



Note that in the above, the x\ and x 2 integrations are over (—6/2,6/2). The above result can be simplified slightly 
by changing variables and writing w — x\ — x 2 and u — x 2 + x\, and noting that for w positive the integration range 
for the variable u is then over (—6 + w,6 — w). Therefore for a generic function f(x, y) even in both its arguments we 
have 

I dxidx 2 f(xi — x 2 , X\ + x 2 ) = I dw / du f(w,u). (39) 
J Jo J-S+w 

Using this change of variables we then obtain 

H b = / dw (6-w) / rdr (r 2 - 2w 2 ) V (yjw 2 +r 2 ^ . (40) 



If we take the limit of small 5 we obtain 

provided the above integral is finite and providing that the integration by parts carried out on the last term of Eq. 
(|35|l is valid. We note that this result is in agreement with Eq. I|29|l as the effective interaction between unit areas 
for thin shells, i.e. for small 6, is 6 2 V. 

The result Eq. I|40|l yields a bending rigidity (here we explicitly include all the integration limits) 



= Y dw (6-w) j rdr (r 2 - 2w 2 ) V (yjw 2 + r 2 ) . (42) 

The same calculation can be carried out for a sphere and one finds that the total energy is given by 

H = 4tt 2 J r\dr x r\dr 2 sm(9)d9 V [\/ (r x - r 2 ) 2 + 2r x r 2 (\ - cos (0))) (43) 

where the r\ and r 2 integrations are over (R — 5/2, R + 8/2). One now makes the change of variables 

r 2 = 2r x r 2 (1 - cos (9)) (44) 

which gives 

H = 7T 2 J rdrd ri dr 2 [(n + r 2 f - (r x - r 2 f\ V (v >i - r 2 ) 2 + r 2 ) , (45) 

where for large R the limits of the r integration can be taken to be (0, oo). Now writing r\ = R + x\ and r 2 — R + x 2 
then switching to the variables u and w as above we find 

r5 r5—w /* oo . 

H = ir 2 dw / du [(2R + u) 2 - w 2 ] / rdr V ( \/ r 2 + w 2 ) . (46) 

JO J-8+w Jo * ' 

Now carrying out the u integral yields 

H = 8tt 2 R 2 j dw(6-w)J^ rdr V (V?' 2 + w 2 ) 

rS 2 4 4 f°° / \ 

+ 7T 2 / dw(8 - w)[-5 2 - -5w - -w 2 } / rdr V ( \/ r 2 + w 2 ) . (47) 
Jo 3 3 3 Jo ^ ' 



8 



7 



Now to obtain the bending energy we must subtract the bulk energy from this result. The area of the flat membrane 
necessary to form this sphere is of radius given by 

= 4»R J + if (48) 
The energy of a flat membrane of this area is then given by 



Hbuik = \ x A S (R) x2n J o dw2 ( S - w ) J q rdr V (Vr 2 + w 2 ^ 

8tt 2 R 2 + / dw{S-w) rdr V[\/r 2 + w 2 ) (49) 

3 J Jo Jo ^ ' 

In calculating Hb = H — Hbulk we see that the term proportional to the surface area cancels and we are left with 

47T 2 f°° 
H b = — J dw w{5 2 - w 2 ) J rdr V {^/r 2 +w 2 } , (50) 

which thus implies 

2n h + Kg = J dww(5 2 -w 2 ) J rdr V \\/ r 2 + w 2 ) (51) 

The result Eq. 1)50(1 is rather significant as it shows that a spherical vesicle has a non-zero bending energy at large 
R when one takes into account the finite thickness of the membrane. When S is small and the corresponding integrals 
turn out to be finite we may ignore the w dependence in the argument of V in Eq. I|51|) to obtain a formula analogous 
to Eq. |@TJh 



jr-fi^ roc 

2n b + Kg = — — J rdr V(r). 



(52) 



EFFECT OF DIFFUSE VAN DER WAALS INTERACTIONS 



Here we consider the problem where the membrane of finite thickness has a different dielectric constant to that of 
the external media/solvent . We take the dielectric constant of the membrane to be e' while the external medium has 
dielectric constant e. This difference in dielectric gives rise to a thermal Casimir effect which is a static van der Waals 
interaction. In general these interactions are not pairwise, but in the diffuse limit, where e — e' is small, the pairwise 
part of the interaction is the dominant one. 

For a system volume C having dielectric constant e' in an external medium of dielectric constant e, the partition 
function for the thermal fluctuations of the zero frequency Matsubara modes of the electrostatic field is given by 



d[4>] exp (S + AS) (53) 



where 



So = -- / dx(V<^ (54) 

and 

A5'=-^— - / dx(V0) 2 . (55) 
2 Jc 

Note that in S the integral is over all space but in AS the integral is only over the volume C containing the media 
of differing dielectric constant to the exterior. As the action in Eq. (|53|l is quadratic the partition function can be 
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written as a functional determinant, however for general geometries this calculation is rather complicated. In the case 
where AS is small one may carry out a cumulant expansion to second order which gives for the free energy F of the 
system 

F-F = -k B T^(AS)o + ^(AS 2 ) ,?j . (56) 

Here Fq is simply the vacuum free energy in the absence of the membrane G. The subscript indicates that the 
expectation is taken with respect to the vacuum measure So and thus 

<0(x)0(y)) o =G (x-y), (57) 

where Go is the Green's function of the free theory with action Sq and the subscript c denotes the connected part. 
Now we note that the first term of Eq. (|56[) is simply a term proportional to the volume of the system and thus does 
not contribute to the bending free energy. The dilute approximation to the bending free energy is thus 

H = -^f(AS\, c , (58) 

In general if one writes 

AS= [ dxi?[0(x)], (59) 
Jc 

H = -^f d*&r<iZ[0(x)TO(y)]>o lC , (60) 
which means that the effective pairwise potential V is given by 

V(pn) = -fc B T(.R[0(x]i?[<KO)])o, c (61) 
In the case of a difference in dielectric constants we thus obtain the potential 

V(x) = -fc B T^^/3 2 ((V0(x)) 2 (V0(O)) 2 > o , c (62) 
Now using the fact that Go(x) = Go(r) where r = |x| we obtain 



then we obtain 



2 , 4 fdC::\' , ., fd 2 G Q 



((V^(x)^(V0(O)) 2 )o, c = ^[^l +2(^7) (63) 



Here the Green's function Go is given by 

- - 

47re^| 



G o(x)=--^ (64) 



and thus we obtain 



y(x)=-A (65) 



|x|' 
where 

A - 

» SfcsT^, (66) 

where A = (e — e')/2 £ ~ (e — e')/( e + e') to leading order in e — e'. This, as is to be expected, recovers the standard 
form of the unretarded van der Waals interaction. 
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We now examine the bending rigidity of a membrane of finite thickness S with this interaction. We regularize the 
interaction by writing 



V(x) 



(x 2 + a 2 ) 3 ! 



(67) 



where a is a short-scale cutoff. Substituting this into Eq. 1421) the r integration is easily performed to yield, after 
some algebra, 



3fci^A 2 f s 

Kb = — / aw (o - w) 

D47T ,/n 



2a 2 



1 



(a 2 + u; 2 ) 2 a 2 + w 2 



which then gives 



3k B TA 2 , fS 2 + a z 



(68) 



(69) 



We see that, as with the eigen-mode expansion method of 8], for small cut-off a the bending rigidity has the form 

(70) 



Kb 



3k B TA 2 8. 



and that the above result is identical to Eq. (0J) up to a rescaling of the microscopic cut-off a. 
The calculation for the sphere is also straight forward to carry out and one obtains 



2Kb = 



k B TA 2 
32ir 



2 rx2 



which along with Eq. I|69(l gives 



k B TA 2 
64tt 



-In 



hi 



(71) 



(72) 



Note that from Eq. (|71() total bending energy for the sphere is positive for all values of 6 /a. We see from Eq. i|69|) 
that the renormalization of Kb due to static van der Waals interactions is rather weak. Although the result depends 
on the short scale cut-off it does so only logarithmically. Physically realistic values of S and a, corresponding to the 
membrane thickness and a typical dipole size or dipole separation gives at most an 0(k B T) renormalization of Kb 

This is much smaller than the experimental values obtained for Kb which tend to be between 3 — 30 k B T 
However Eq. 1)72(1 predicts a K g which depends strongly on S and it is conceivable that van der Waals interactions 
make a significant contribution to K g . Unfortunately few experimental measurements or estimates exist for K g . 



VI. CONCLUSIONS 



In this paper we have revisited the problem of the renormalization of the bending and Gaussian rigidities of 
membranes by long-range interactions. These renormalizations may be calculated via a geometric approach applied 
to cylindrical and spherical geometries. The result obtained for Kb is found to agree with that found for a perturbative 
analysis about a flat membrane in a general approach proposed by Netz 0. We also rederived the Netz result for 
a fiat membrane by considering an ensemble where the membrane is thought of as being drawn from a reservoir of 
flat non-fluctuating membrane showing the equivalence of the two approaches from an ensemble point of view. Using 
the geometric approach we obtained the general result for an infinitesimally thin membrane that K g = —2Kb, the 
somewhat surprising result that the Gaussian rigidity is renormalized by a long-range potential with a sign opposite to 
the bending rigidity. This effect has also been seen in more specific mean-field studies of charged membranes. We then 
derived analogous formulas for the bending and Gaussian rigidity when the membrane has a finite (but small relative 
to the radii of curvature) thickness. Finally we calculated the bending and Gaussian rigidities induced by the thermal 
Casimir force, or static van der Waals interactions, for a dilute system. In agreement with our previous studies we find 
a negative contribution to the bending rigidity but with exactly the same functional dependence on the thickness 5 of 
the membrane and the microscopic cut-off a when a is small compared to S. This agreement is reassuring as although 
microscopic details are dominating the physics they are doing so in a rather universal way which is insensitive to the 
regularization scheme being employed. The Gaussian rigidity is found to undergo a positive renormalization due to 
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van der Waals interaction. This renormalization strongly divergent as a — > behaving as S 2 /a 2 . As pointed out in g|, 
van der Waals interactions only weakly favor the formation of tube like structures, such as t-tubules and it is unlikely 
that they can stabilize cylindrical geometries thermodynamically. However if these structures are formed via another 
physical or biological mechanism, then attractive van der Waals forces may contribute to their stability in that they 
will impede the formation of spherical budding instabilities and thus enhance the metastability of these structures. 
Whether this enhancement of metastability is significant depends the value of 8/a and contributions of a similar 
functional form which will come from the non-zero frequency frequency Matsubara modes which are responsible for 
the fluctuating part of the van der Waals interaction. 

A final comment on the conflicting results Eq. JSJ) and Eq. © for the renormalization of the bending rigidity of 
a salty membrane is in order. In the case of a system with salt outside, the bare interaction induced between the 
membrane components, in the diffuse limit where the salt concentration within the membrane is small, behaves as 
exp(— 2mr)/r 2 where m is the external Debye mass. In this limit the two results should give the same result as the 
interaction is sufficiently short range. In the Netz approach an effective pairwise interaction is computed, however in 
the eigen-mode expansion of Pincus and Lau it is clear from their calculation that they are summing all terms in the 
eigen-mode expansion and not just those which correspond to the dilute-pairwisc limit. Concretely the logarithms in 
the expressions given in |6( are not expanded to second order. This means that their result is inherently taking into 
account multiple scattering events and thus non-pairwise interactions. All the same one should bare in mind that the 
physical situation where one has salt in the membrane but non outside is rather unlikely. 
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